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In the standard hot big bang nucleosynthesis (BBN) mo delB 0, the pri- 
mordial abundances of 1 H, 2 H, 3 He, 4 He, and 7 Li fix the baryon density of the 
universe, Q^, via the baryon-to-photon ratio, 77, for a given Hubble parameter. 
Recent observations of Li show@ that its intrinsic dispersion in metal-poor stars 
is essentially zero, and the random error in the mean Li abundance is negligible. 
However, a decreasing trend in the Li abundance towards lower metallicity, plus 
6 Li detections^ indicate that its primordial abundance can be inferred only 
after allowing for nucleosynthesis processes in the Galaxy (Galactic chemical 
evolution, hereafter GCE). We show that the Li vs Fe trend provides a tough 
discriminant between alternative models for GCE of light-elements. We criti- 
cally assess current systematic uncertainties, and determine the primordial Li 
abundance within new, much tighter limits. We show that the Li constraint 
on Qb is now limited principally by uncertainties in the nuclear cross-sections 
used in BBN calculations, not by the abundance itself. A clearer understanding 
of systematics allows for a much more accurate inference of the primordial Li 
abundance, sharpening the comparison with 4 He and deuterium and the result- 
ing test of BBN. We find that the Li data are in good agreement with 4 He and 
with "high" deuterium values, but that low deuterium abundances are at best 
marginally within the Li range. 

An important test of BBN is concordance between the observationally-inferred primordial 
abundances of the light elements. However, all current values involve considerable systematic 
uncertainties. Estimates of the 4 He primordial mass fraction, Y p , had settled® around 
Y p = 0.230 ± 0.005, but different systematics and underlying stellar He I absorption may 
imply a higher valueO near 0.245. For deuterium, quasar absorption line measurements give 
both "low" SB abundances around D/H = (3-5) xl0~ 5 and "high" 00 values around 
D/H = (15-25) xl0~ 5 . It is unclear which value (if either) represents the correct primordial 
value. 3 He presents even greater difficulties associated with its highly uncertain yield in 
low-mass stars 0' 0, and currently provides unreliable constraints. Recent observations!! 
of Li have greatly improved, therefore we examine the random and systematic uncertainties 
associated with the primordial abundance and its interpretation with respect to BBN. 

In inferring the primordial 7 Li abundance, A(Li) p (where A(Li) = logio(n(Li)/n(H)) 
+ 12.00), from observations of metal-poor stars, systematic errors arise in: (1) the assessment 
of Li GCE prior to a given star forming; (2) the correction for depletion of a star's initial 
surface Li; (3) the measurement of the current abundances; and (4) possible confusion by 
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anomalous objects. We examine each factor below, and summarise its impact in Table 1. 

GCE: The GCE contribution to 7 Li was long believed to be negligible for metal-poor 
stars. However, observations of very metal-poor stars showing (1) lower Li abundances^ 
and (2) measurable 6 Lifl (which comes solely from Galactic cosmic ray (GCR) reactions), 
show that GCE cannot be ignored. Li GCE can be constrained either empirically from the 
observed evolution, or by modeling the sources and sinks of the element. Here we pursue 
both routes, partly to investigate the range of possible values, and partly to learn about 
GCE itself, since an accurate Li trend can constrain not only standard GCR and //-process 
nucleosynthesis, but also other mechanisms which may produce Li in Population II, e.g., 
those suggested to produce primary beryllium along with Li. Sources of lithium in the 
oldest (Population II) stars are GCR nucleosynthesis of 6 Li and 7 Li, and the supernova v- 
process which produces 7 Li and n B. The sink is astration — the destruction of Li in stellar 
interiors subsequently ejected into the interstellar medium (ISM). 

Previously^ we examined GCE empirically, tracing the production of Li as the iron 
abundance increased, as a regression in logarithmic abundances: 

A(U) = a + /?[Fe/H] (1) 

where [Fe/H] = log 10 (Fe/H) star — log 10 (Fe/H) sun . (Iron is a convenient and widely-measured 
diagnostic of GCE. It is always produced and never destroyed by stars, and thus serves as a 
chronometer.) We obtained values in the range (3 = 0.07-0.16, depending on the adopted 
[Fe/H] values, the errors, and whether some points are excluded from the fit. 

In the present work, we also investigate a fitting form which better follows Li production 
by GCE, and which simplifies the extrapolation to the primordial value. Li production is 
proportional to the cumulative number, Nsn, of Type II supernovae, as these are both GCR 
accelerators and the site of the //-process. It thus is important to establish the primary tracer 
of such supernovae. If the cumulative supernova rate is well reflected by the iron abundance 
(Nsn oc Fe), then a fit to linear abundance scales is appropriate: 

Li/H = a' + b'Fe/Fe Q (2) 

Here a' directly measures the primordial 7 Li abundance (in the absence of systematic errors), 
while b' probes GCE. The linear fit parameters are sensitive to systematic Li errors; a change 
by A cal dex in the log shifts both a' and b' by a factor 10 Acal . We find a' = 1.0-1.2 xl0~ 10 
and b' = 40-180 xl0~ 10 . 
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On the other hand, if oxygen (which is more difficult to measure than iron) is a better 
tracer of Type II supernovae than iron, then N$n oc O, and we expect: 



Li/H = a + 6O/O (3) 

where O/O = (Fe/Fe Q ) 1+a; . Recent observations liil li^C show uo = —0.31. In this case the 
data indicate a = 0.9-1.2 xl0~ 10 and b = 9-34 xlO" 10 . 



We also compute the Li-Fe trend expected from a one-zone (closed box) GCE mo deio m @ 
which includes GCR and stellar nucleosynthesis, to compare with the observational results. 
(See Figure 1 caption for details.) Figure 1 shows the different Li components for the model 
with 7 Li p = 1.23 x 10~ 10 . We fit the model by regressions of the forms (P - @, over 
the metallicity range of the recent dataB, and find b = 8.6 x 10~ 10 , V = 4 x 10~ 9 , and 
P = 0.03 - 0.07 (for input 7 Li p = (0.9-1.9) xl0~ 10 ). (The model's "linear slopes", b and 
b', are independent of the input 7 Li p , while its "log slope", /3, does depend on 7 Li p .) This 
model's Li GCE slopes (6, V and (3) are at the lower range of those found for the data, 
and for a closed-box model with only GCR and the //-process. The poor agreement can 
be alleviated somewhat by manipulation of free parameters. For example, models including 
outflow of supernova ejecta (open box) and/or the inclusion of Li yields from AGB stars 
with M < 5M produced steeper slopes, b = (11 - 17) x 10~ 10 , b' = (6 - 9) x 10" 9 , and 
(3 = 0.06 — 0.08 (for various types of models with 7 Li p = 1.23xl0 -10 ), closer to the range 
seen in the data. This demonstrates the power of a reliable observational Li versus Fe trend 
to constrain the form and parameters of GCE models. 

While comparison of model slopes with the observations can teach much about GCE, to 
infer the primordial Li we use the observed slopes (a procedure very similar to that used 
for 4 He). From the linear fits to the data and our previous analysisB, we estimate that the 
GCE contribution to this metal-poor turnoff sample is —0.11^09 in the log (see Table 1). 

Stellar Depletion: Stars burn Li, preserving at most a thin outer shell containing a few 
percent of the star's mass. Possible in situ depletion of Li has long been regarded the major 
systematic uncertainty in inferring A(Li) p from present-day abundances. Stellar evolution 
models predict the depletion factors. The simplest models imply almost no destruction 
(<0.05 dex, possibly <0.01 dex) in very metal-poor turnoff stars^. Models incorporating 
rotationally-induced mixing had predicted large depletion factors ~ 1 dex, though more 
recent efforts give lower values" ~ 0.2-0.4 dex, and predict a range of depletion factors 
from star to star. The negligible intrinsic spread found for very metal-poor turnoff stars, 
< 0.02 de X , ra ,e S out rotational depletion even a S low aa 0.1 dexi. A s d lffuS1 on is al S o 
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absented, we conclude that in situ depletion is minor, <0.1 dex, and possibly as little as 
~ 0.01 dex. 

Abundance Analyses: A Li abundance is derived via a parameter- and model-dependent 
analysis of a stellar spectrum, and systematic uncertainties propagate through to A(Li) p . 
Effective-temperature calibrations can differ by up to 150-200 K, higher temperatures re- 
sulting in higher Li abundances by 0.065 dex per 100 K. The scale initially adopte di gives 
temperatures cooler than a more recent calibration @ by on average 120 K. We now ad- 
just the abundances (Table 1) to the newer calibration^ , but note that systematic errors 
of ±120 K may still exist. This is one of the largest contributions to the uncertainty in 
A(Li) p . Fortunately, errors in the surface gravity, microturbulence, or damping parameters 
are neglig iblei. Concerns about 1-D, plane-parallel model atmospheres have been reduced 
by simulations of solar-type granulation® which show that the Li abundance is underes- 
timated in the 1-D approximation by < 0.10 dex, and possibly < 0.01 dex, depending on 
the theoretical prescription for microturbulence. Consistent results in the metal-poor star 
HD 140283c3 from the Li 6104 A and 6707 A lines inspire further confidence. However, 
models with greater convective flux can lead to Li abundances higher by 0.08 dexED. Cor- 
rections for non- LTE© are only —0.01 to —0 03 dex, and the uncertainty in the ^/-values 
is only 0.02 dex (ltr)Hl. 

Anomalous objects: Apart from the grossly Li-depleted star G186-26, only one of the 
remaining 22 objects in our sample was rejected by outlier-detection algorithms, changing 
the mean abundance by only ~ 0.005 dex. Similar un-recognised objects would affect the 
result by < 0.01 dex. 

We can use the primordial element abundances to fix the one free parameter of the 
standard hot Big Bang nucleosynthesis (BBN) mo deli I, the baryon-to-photon ratio, r). 
From this, the baryon density of universe, Qb, may be deduced (for a given Hubble pa- 
rameter, h). The inferred primordial abundance for Li (Table 1) is A(Li) p = 2.09lo!i3 
(Li/H = 1.23lo;32 x 10~ 10 ), where the errors incorporate statistical (negligible) and sys- 
tematic (more significant) effects. These errors are now sufficiently small that the range of 
corresponding rj values (see below) is dominated by the uncertainties in the input nuclear 
cross-sections used BBN calculations rather than in the abundance. (Uncertainties in BBN 
give rise to a range of r\ values similar to those quoted below even for a perfectly determined 
value of A(Li)p.) 

An important test of the BBN model is whether the inferred primordial abundances give 
concordant values of rj. This is best tested by establishing likelihood distributions (as a func- 
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tion of rj) for each element, convolving the theoretical and observational uncertainties 00. 
Figure 2 shows the likelihood distributions for 4 He and four possible values of the primordial 
7 Li abundance, all of which give excellent agreement. Overall concordance (at 95% CL) 
occurs for r] 10 = 10 10 r] = (1.4-4.9), (1.5-4.4), (1.7-3.9), and (1.8-3.6), for 10 10 x 7 Li/H = 1.9, 
1.6, 1.23, and 0.9 respectively. We can then use this result to assess the diverse deuterium 
values. For high D/H (2.0xl0~ 4 ), the peak of the D/H likelihood function (not shown) is at 
7/10 = 1-7, (95% CL = 1.4-3.8), in very good agreement with the results from 4 He and 7 Li. 
For low D/H (3.4xl(T 5 ), the peak of the D/H likelihood function is at r] 10 = 5.2 (95% CL 
= 4.6-6.1), which would require 7 Li p at the upper end of the possible range. However, if the 
low D/H was even only slightly higher, at SxlO -5 ^, then the D/H peak occurs at r] 10 = 
4.0 (95% CL = 3.6-4.6), consistent with the ranges for 4 He and 7 Li. 

The baryon density corresponding to r\ = (1.7-3.9) xlO -10 is Qb =(0. 025-0. 057)//i 2 
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(where is the Hubble constant in units of 50 km s 1 Mpc 1 ), significantly below the 
lower 95% CL on Q m = 0.251^^, thus maintaining the requirement for non-baryonic 
dark matter. 
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Table 1 Inferred Primordial Lithium Abundance 



Observed:§ (A(Li))_ 2 . 8 = +2.12 ±0.02 

Corrections to apply (logarithmic): (Estimated 

(1) GCE/GCR: Uncertainty) 
previous analyses @ —0.14 to —0.05 

log data fit (eq. (1)) -0.20 to -0.09 

linear data fit (eq. (2)) -0.12 to -0.04 

linear data fit (eq. (3)) —0.16 to —0.05 

model fits (eq. (2)-(3)) -0.12 to -0.02 

Adopted (excludes model): —0.11 

(2) Stellar depletion +0.02 
(3a) T eff scale zeropoint +0.08 ±0.08 
(3b) 1-D atmosphere models +0.00 itoloo 
(3c) Convective treatment +0.00 ±g;gg 
(3d) NLTE -0.02 ±0.01 
(3e) ^/-values +0.00 ±0.04 
(4) Anomalous objects +0.00 ±0.01 



Total -0.03 
Inferred: A(Li) p = +2.09 



+0.07 
-0.09 
+0.08 
-0.02 



-0.19 
-0.13 



+0.19 
-0.13 



Table Caption: 

Table 1: Weighted mean and 95% CL uncertainty of observed Li abundances for a very 
metal-poor turnoff sampleS with ([Fe/H]) = —2.8, and corrections required to deduce the 
primordial value. (The weighted mean differs slightly from the robust mean, viz. 2.110.) Five 
estimates of the logarithmic correction for GCE are listed, based on the previous analysisS 
(logarithmic fits and observed 6 Li/ 7 Li ratio) and the new work in this paper (linear fits 
and various GCE models). The model fits are based on the model slopes and the inferred 
deviation of the primordial value to the weighted mean at [Fe/H] = -2.8. The various error 
estimates, which include random and systematic uncertainties, are clearly non-Gaussian, so 
combining them is an imprecise and subjective process. We take quadratic sums for the 
positive and negative uncertainties separately, and regard these as estimates of the 95% 
confidence limits. 
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Figure captions 

Fig. 1 Contributions to the total predicted lithium abundance from the adopted GCE 
modelEU , compared with low metallicityl^ and high metallicity® stars. The solid curve 
is the sum of all components. Additional stellar production mechanisms of 7 Li are required 
for stars with [Fe/H] near zero, but are expected to be unimportant for the lowest metallicity 
objects. The primordial component of 7 Li decreases at high metallicity due to astration, but 
other components increase with metallicity as described in the text. The GCE model 03 EB 
is a one-zone (closed box) model which includes GCR and stellar nucleosynthesis. The 
latter does not reproduce the observed O-Fe scalingEH , so we retain the calculated O 
evolution and force Fe to match [O/Fe] = a; [Fe/H], (with uj = —0.31). For GCR production 
of Li, the model assumes that the cosmic ray flux is proportional to the supernova rate, and 
that GCR abundances match the ISM. These assumptions fix the linearity of Li-0 scaling 
at low metallicity where a + a dominates. The scale factor is set by the GCR particle 
spectrum and confinement, for which we take a source spectrum oc p~ 2 , and an escape 
path-length A = 100 gem -2 . The model requires three remaining inputs: (1) an adopted 
primordial 7 Li p abundance, (2) the overall normalization for all GCR production, and (3) 
the //-process contribution. As GCR nucleosynthesis also produces beryllium and boron, 
we use the meteoritic 9 Be and 10 B abundances to establish to overall normalization, which 
then fixes the GCR contributions to 6 Li and 7 Li. With these normalizations, the modeled 
evolution of Be and B fits available Population II observations 0' EDI. 

Fig. 2 Likelihood distributions for four values of primordial 7 Li/H (10 10 x 7 Li/H = 1.9 
(dashed), 1.6 (dotted), 1.23 (solid), and 0.9 (dash- dotted)), and for 4 He (shaded) for which we 
adopt Y p = 0.238 ± 0.002 ± 0.005 (random and systematic uncertainties) 01 For 7 Li/H = 
1.6xl0~ 10 (v4(Li) = 2.20), there are two likely values of 7710 = 10 10 i] = 1.9 and 3.6, because 
the predictions are not monotonic in 7 Li. For 7 Li/H < 1.1 x 10~ 10 (A(Li) < 2.04), the 
Li abundance is at or below the BBN predicted value, so there is only one peak, at r] w 
~ 2.6; uncertainties in the prediction and observation prevent the likelihood function from 
vanishing. The peaks of the combined distribution (the product of Li-R^rj) and Lruiv)] not 
shown) are at roughly the same value of r] as in the individual L-7 Li (r\) distributions. 
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Figure 1 
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